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Multifractal  stationary  random 

processes 


Random  cascades  (Mandelbrot) 
Intuitive  construction 


Theoretical  construction 


Some  particular  examples 


Random  cascades  (Mandelbrot) 


Dyadic  case 

•  ln  =  TXn,  A  =  1/2 

.  /  =  [0,  T[,  I0  =  [0,  T/2[,  /oi  =  [T/4,  T/2[,  . . . 

•  Mln(ISo...Sn)  =  Wg0,,,SrlMlrl_1  {Iso...8n-l  ) 

•  (WS0...Si }  are  i.i.d, 

•  E  (WS0...Si)  =  1/2  (=>  Min  martingale). 

•  “Discrete  scale  invariance”  (around  0) 

Jo /2  Min (dt)  =l™  W  &  Mln_±{dt) 

Mln(dt/ 2)  =Io“  WMin_1(dt) 

Mln  {dtln  /  ln—\)  ="  WMin_1  (dt) 

By  taking  n  — ■>  -foo 

rt/2  7  rf 

/  M(dt)  =  W  /  M{dt) 

J  o  Jo 


>  NO  :  not  stationary,  no  continuous  scale  invariance 


Towards  “continuous”  cascades 


Discrete  cascades  (dyadic  case) 


Mtjdtln/ln-t)  l^WMln^(dt) 


=>  Try  to  represent  it  as  a  discretization  of  an  underlying 
continuous  construction. 

•  Mv{Vdt )  =law  Wi>,i(t)Mi(ldt),  V  <  l  <T 

•  law  of  Wy:i(t)  depends  only  on  V jl  (and  i.i.d). 
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Why  infinitely  divisible  laws? 


•  Mi(dt)  =law  Wi,i'(t)Mi'{dt ),  l  <  V  <T 

=>-  In  Wi'  jp  ~  In  Wi^i  +  In 

=>  In  Wi^t  =  sum  of  2  i.i.d  variables  (for  l  =  \JVT  ) 
ui(t)  =  In  arbitrary  number  of  i.i.d  var. 

=>  u>i(t)  is  infinitely  divisible. 


•  P(dt,  dl )  stochastic  “infinitely  divisible  noise”  on  the 
half-plane  (t,  Z)  G  1 R  x  1R+*  with  respect  to  the 
measure  fi(dt,dl), 

-Ac  <S+,  P(A)  infinitely  divisible 

-  the  law  of  P(A)  depends  only  on  fi{A) 

-  P(A)  and  P(B)  independent  iff  An  B  =  0 


2d-representation  of  discrete 

cascades 


^ In  (s0  •  •  •  sn)  —  X/i=0  HWso...ai) 


1 


M1 

*2 

*3 


Time  t 


(Jtri  (s0  . . .  sn)  =  P(A2-n  (s0  . . .  sn)) 
lnWSo...Si  i.d  =^>  fi(dt,dl)  =  dtdl/l 2 
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Continuous  cascades 


Q) 

(0 

O 

CO 


T/2 

T/4 

T/8 


A1/8(000) 


a1/9(010) 


InW, 


\^nW0 


n^V„ 


M, 


M, 


M. 


Time  t 


=  \nWhT{t)  =  P(Mt) 


Mi (dt)  d=  e“l{t)dt 

>  Mi>  (dt)  =law  Wi',i(t)Mi(dt),  V  <  l  <T 
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Infinitely  divisible  laws 


•  Definition 

X  is  infinitely  divisible  iff  it  is  a  sum  of  an  arbitrary 
numbers  of  i.i.d.  variables 

•  Examples 

—  Gaussian  variable  (mean  m  ,  variance  a2  ) 


p(x) 


4>x(q)  =  InE  ( eiqX^j  =  imq  —  q2(J2 / 2 
<fx  (q)  is  the  cumulant  generating  function  of  X 
—  Poisson  variable  (intensity  A  ) 
p(x)  =  e~xXk /k\ 


fx(q)  =  X(etq  -  1) 


•  Theorem 


A  sum  of  independent  infinitely  divisible  variables  is  an 
infinitely  divisible  variable. 


Moreover  faX+pY^q)  =  (fx(aq)  +  <f>y(l3q). 
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The  Levy-Khintchine  formula 


•  Combining  gaussian  and  (an  infinite  number  of)  poisson 
variables  : 

0(g)  =  imq  —  q2a2 /2  +  f  ( eiqx  —  1  )v(dx) 
u(dx)  ~  “intensity  measure” 
a  priori,  total  intensity  bounded:  f  P(dx)  <  +oo 
actually,  intensity  around  0  can  be  +oo  ( x~2dx ) 


•  Levy-Khintchine  representation 

f  eiqx  —  1  —  iq  sin(x) 

0(g)  =  imq  +  /  - — - 

— »  v(dx)  is  the  Levy  measure 
->■  fly  v{dx),  f_fo  v(dx) / x2 ,  fy°°  v(dx)/x 2  <  +oo 

—  v(dx)  =  5{x)  — >  gaussian 

—  v(dx)  =  S(x  —  xo)  -^poisson 

—  J  is(dx)/x2  <  +oo  —>  compound  poisson 

—  a-stable,  gamma,  Student,  . . . 
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“Infinitely  divisible  noise” 


Definition  ( Rosinski ;  1989 ) 

P(dt,  dl )  is  an  independently  scattered  infinitely  divisible 
random  measure  distributed  on  the  half  plane  <S+  with 
respect  to  the  (deterministic  measure)  fi(dt,  dl) 

•  {An}n  disjoint  sets  of  S+  ,  {P(An)}n  are  i.i.d. 


oc 


P  (ur=1An)  =  p(An)  ,  a.s. 


n—  1 


•  P(A)  is  infinitely  divisible  with 


Notation  :  if(q)  =  <f(—iq)  whenever  it  is  possible  and 
f>(q)  =  +oo  otherwise  (f>(q)  is  convex) 


f>(q)  =  +00,  if  E  (eqP(A))  =  +oo, 
E  {eqP<K^)  =  e^(qAAf  otherwise. 
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Continuous  cascades 


Mi(dt)  =def  ep(-Al^dt 

•  Schmitt  F.  and  Marsan  D.,  2001 

-  basic  ideas  for  construction  of  ui(t)  =  P(Ai(t)) 

-  scaling  of  e UJl  versus  l 

•  Bacry  E.,  Delour  J.  and  Muzy  J.F.,  2001 

-  Id  representation  of  the  log-normal  case  (MRW) 

-  exact  power  law  scaling  of  the  moments 

-  no  proof  of  convergence 

•  Barral  J.  and  Mandelbrot  B.B.,  2002 

-  full  2d  compound  poisson  construction 

-  non  degeneracy,  Lp  convergence,  . . . 

-  multifractal  formalism 

•  Bacry  E.,  Muzy  J.F.,  2003 

-  full  2d  log-inf.  divisible  construction  (MRM,  MRW) 

-  non  degeneracy,  Lp  convergence,  . . . 

-  exact/asymptotic  scaling  of  the  moments 

•  Abry  P.  ,  Chainais  P.,  Riedi  R.  2003 

-  full  2d  log-infinitely  divisible  construction  (IDC) 

-  L2  convergence 

-  non  power  law  scaling  of  the  moments 


li 


Multifractal  Random  Measures 

(MRM) 

(Bacry,  Muzy,  2003) 


•  Ai(t)  domain  : 

-  General  definition 

=  v>i,  \t'-t\<f(o/2}. 

fi°°  f(s)/s2ds  <  oo  and  /(/)  —  l  for  l  <  L 

-  Definition  for  “exact  scaling” 

=  l  for  l  <  T  and  =  T  for  l  >  T 

•  ui(t)  =  P(Ai(t)) 

-  Levy  measure  v{dx) 

-  ^(1)  =  0  {eB1^  martingale) 

-  3e  >  0,  i)(l  +  e)  <  +oo  (E  (gU+eVdL)  <  -boo) 

•  Mi(dt)  =  e^l^dt 

•  M(dt )  =  lim/_^0+  M(dt ) 
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Existence  of  the  limit  MRM  measure 

(Bacry,  Muzy,  2003) 


Theorem  1  There  exists  a  measure  M(dt)  such  that 
(i)  with  probability  one ,  Mi(dt )  — >  M(dt)  weakly 
(ii)  Mt  g  JR,  =  0, 

(in)  for  any  bounded  set  K  of  JR, 

-  M(K)  <  +oo 

-  E  (M(K))  <  \K\. 

Proof 

•  =  0  E  (eP1^}  =  1. 

•  is  a  left  continuous  positive  martingale 

•  Kahane  J.P.,  Chi.  Ann.  of  Math.  8B;  1-12  (1987) 
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Continuous  scale  invariance 

(Bacry,  Muzy,  2003) 


Mv{l'dt)  =  Wv,i{t)Mi{ldt),  V  <  l  <T. 

=>  M\i(Xldt)  =  W\(t)Mi(ldt),  V  =  XI, 
=>  M\i(Xdt)  =  W\(t)Mi(dt), 

=>  Mxl([0,Xt])  =  Wx(0)Ml([0,t]). 


On  the  one  hand 
MXi([0,Xt}) 


ew*,(uW  =  x 


f  ew^Xu)du 
0 


On  the  other  hand 


Wx(0)Ml([0,t})  =  Wx(0) 


ew,(uW 


Consequently  we  want 

Aeu.Al(A t)t™Wx  f  eMu)du 

J  0 

Or, 

toxi(Xt)  l^’ln(Wx/X)+col(t)m 
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Scaling  of  uji(t) 

(Bacry,  Muzy,  2003) 

We  want  {uj\i{\t)}t<T  =law  ln(WA/A)  +  {ui(t)}t<T, 

Lemma  1  (Characteristic  function  of  uoi(t)) 

]j^  j  —  —  l  — l  PlOk~tj) 

where 

-  pi(t )  =  /i(*4/ (0)  n^(t).), 

-  E?=iEi=i«(j'^)  =  ¥>(ELi^)- 

□ 

The  exact  scaling  domain  f^e\l)  is  the  only  domain  which 
satisfies 

P\i( At)  =  —  log  A  +  p/(t),  l  <  T,  A  <  1,  t  <  T 

In  this  case,  ln(ITA/A)  is  infinitely  divisible  (indep.  of  M) 
with 

£  /^e^ln(WA/A)\  _  ^-(p(q)' 


15 


Continuous  scale  invariance  of  M^e\t) 

(Bacry,  Muzy,  2003) 


Theorem  2  (Continuous  invariance  of 

{M^([0,M})}t<T  la^Wx{M^([0,t})}t<T, 
where  ln(W\/X)  is  infinitely  divisible  (indep.of  M )  with 


Theorem  3  (Moment  scaling  of  M^e\t)) 

E  (M(e>([0,«])9)  =  (2)  E  (M(e)[0,r])9)  ,  Vi  <  T. 
where 

C  (q)  =  q-,>P(q) 
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Moments,  Degeneracy  of  M^e\t) 

( Bacry ,  Muzy,  2003) 


Theorem  4  (Existence  of  the  moments  of  M(e\t)) 
Let  q  >  0  then 

(i)  ((q)  >  1  ==>  E  (M^  ([0,  £])9)  <  H-oo  and 
sup^  E  (m^  ([0,  t])q^j  <  +oo. 

(ii)  if  7^  0,  E  (M^e^([0,  t])q)  <  +oo  =>  ((q)  >  1. 


Proof  for  (n) 

E  (M^([0,t])A  =  E  ((M^([0,t/2])  +  M<ye\[t/2,t}))q\ 

>  E  (M^e\[0,  t/2])q\  +E  (M(e){[t/2,t})q'\ 

Using  scaling  of  the  moments 

E  (M(e)([0,t])d  >  21_C(9)E  (^M(e)([0,t])h  , 

and  consequently  ((q)  >  1. 

□ 


Theorem  5  (Non  degeneracy  of  M^e\t)) 

(H)  3  e  >  0,  £(1  +  e)  >  1 

if  (H)  holds  then  E  (M^e^([0,  t]))  =  t. 
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Summing  up  results  on  M^e\t) 

(Bacry,  Muzy,  2003) 


(i)  ^(1)  =  0  =>  existence  of  M ^ 

(ii)  =3  e  >  0,  C(1  +  e)  >  1  =>  E  (M(e)  ([0,  £]))  =  t 

(iii)  (( q )  >  1  “  ”  E  (M^e^  ([0,  £]))  <  +oo 

(iv)  {M^([0,Xt])}t<T  =lawWx{M^([0,t})}t<T, 

(v)  E(M<e)([0,t])«)  =  (±)aq)E(M^[0,T])‘i)  ,  Vi  <  T. 
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From  M^it)  to  M(t) 

(Bacry,  Muzy ,  2003) 


Theorem  6  (Degeneracy,  asymptotic  scaling  and 
moments  of  positive  orders  of  M(dt )) 

(i)  M^e\dt)  =a-s-  0  4=^  M(dt )  =a-s-  0, 

Moreover ,  if  M(dt)  ^  0,  one  has 

(ii)  with  probability  one , 

Vt  >  0,  M([0,  At])  -  At])  ;  when  A  0+. 

(Hi)  E  {M^e\[ 0,  t])q )  <  +oo  <=>  E  (M([ 0,  £])g)  <  +oo. 

fW  E(M([0,i])«)~t^o+  (^)C(9)E(M([0,T])«). 

□ 
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Multifractal  Random  Processes 

(MRW) 

(Bacry,  Muzy,  2003) 


•  Definition  (subordinated  process) 


X^(t)  =  B(M([0,t}))- 


Equivalent  definition  (stochastic  integral) 


where 


X(t)  =  lim+  Xi(t), 


Xi(t)  =  [  eul(u^2dW{u) 

Jo 
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Main  theorem  on  MRW 

(Bacry,  Muzy ,  2003) 

Theorem  7  Under  hypothesis  £(1  +  e)  >  1  (non 
degeneracy  of  M ), 

(i)  C (q)  >  1  =>  E  ( \X(t)\2q )  <  +oo. 

(ii)  E  (|X(t)|2g)  <  +oo  ==>  ((q)  >  1. 

(Hi)  {X^(t)}t<T  =lawWx{X(e\t)}t<T, 

(iv)  E  (\X(e\t)\2q)  =  (±)C{q]E  (\X(e\T)\2q)  ,  Vt<T. 

(v)  E  (\X (t)\2q)  (^)C(9)E(|X(T)|2g)  . 
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Numerical  simulation  of  MRW 

( Bacry ,  Muzy ,  2003) 


•{e[k]}k£7L  Gaussian  white  noise,  variance  1. 

•  In  =  2~n 

•  Xin(t)  =  J2ktoa'/Le  !,i<2  ~e[k\. 

Theorem  8  (Convergence)  If  (( 2  +  e)  >  1  then 

lim  {Xln{t)}tl={X(t)}t. 

n— >+oo 

Proof 

liniji^+oo  Min  (dt)  =m-s ■  M(dt) 
where  Min([0,t})  = 

□ 


=>  X(t)  is  the  continuous  limit  of  a  random  walk  with 
stochastic  variance 
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Log-normal  MRW 

( Bacry ,  Muzy,  2003) 

•  Tz(dx)  =  2X2S(x), 

•  ^(q)  =  (lm  +  A2g2  with  m  =  —A2. 

•  C(q)  =  di1  +  a2)  -  AV- 


a2  —  1,  T  =  512  and  A2  =  0.025. 


•  Id  construction  :  Bacry ,  Delour ,  Muzy ,  2001 

^  ^ Z(t)  /  \ 

Xi(t)  =  Xi(t  —  l)  +  ayle^~ ei(t ),  (step-wise  process) 
E(uji)  =  —Var(uJi)/2 


Cov(uii  (0), 


4A2ln(r^iy)for  \t\  <T 
0  otherwise 


X(t)=]iml^0+Xi(t) 
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Log-Poisson  MRW 

(Bacry,  Muzy,  2003  and  zz  Barral,  Mandelbrot,  2002) 


•  v(dx)  =  ^5{x  —  xo),  with  7  =  Axq,  =  In  5 

•  ((q)  =  qrri  +  A(1  —  5q)  (m  such  that  £(1)  =  1) 


cr2  =  1,  T  =  512,  A  =  4  and  cc0  =  0.05. 

•  Log-Poisson  compound 

f  is(dx)x~2  —  C  <  -f-oo,  pin  w  =  v(dx)x~2 / C 

Compounded  with  In  VP,  intensity  C 

C(q)  =  qTn  +  (7(1  —  E  (VP9)),  (m  such  that  £(1)  =  1) 
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Other  MRW’s 

( Bacry ,  Muzy ,  2003) 


•  log-ce-stable  MRW 

left-sided  ct-stable  density  : 


I  C\x\x  a  if  x  <  0 
v{dx)  =  < 

I  0  if  x  >  0 

with  C  >  0  and  0  <  a  <  2 
C(q)  =  qm  —  cra|g|a 
stable  =>  Id  construction 


•  log-gamma  MRW 

v(dx)  =  Cry2xe~lxdx  for  x  >  0 
£(<?)  =  qm-  C72  In  with  7  >  1 


•  log-Stndent,  log-Pareto,  . . . 
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Log-normal  MRW 

( Bacry ,  Muzy ,  2001,  2003) 


•  Continuous  limit  of  a  rand,  walk  with  stoch.  variance 
V(t)  =  —  l)  +  a2vle^~ ei(t),  (step-wise) 


E  (oji)  =  -Var(u)i)/2 
Cov(uJi(0),un(t))  =  | 

X(t)=  lim^0+  Xi{t) 
•  C (q)  =  q-q(q-  i)A2 


4A2ln(rFtw)for  1*1  -T 

0  otherwise 


•  Only  3  parameters  : 

-  a2  :  deterministic  random  walk  variance 

-  T  :  correlation  length 

-  A2  :  “multifractal  coefficient” 

A  =  0  =>X(t)  =  B(t) 
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log2M(q,l)  MRW  sample 


Multifractal  analysis  of  a  MRW 

( Bacry ,  Delour,  Muzy ,  2001) 

800 
600 
400 
200 

1990  1995  2000 

0 

-10 

-20 

5  10 

log2(l) 

2 

0 

-2 


time  (year) 
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In  P^Sx) 


Moments  and  law  deformation 

across  scales 

( Bacry ,  Muzy,  2001,  2003) 

•  Explicit  g-order  moment  formula 

E  (X(t)‘ 


•  q  >  1/A2  ” 


E(X(t)2<?)  =  +oo 


-10 


-5 


15 


MRW 


Numerical  estimation 


log-normal  approx. 


20 


-10 


-5 


10 
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Discretization 

(Bacry,  Kozhemyak ,  Muzy,  2004) 


Modeling  discrete  data  with  stationary  increments 

-  Discrete  model  Xa  [n\  =  X(nA) 

-  Stationary  incr.  5Xa[ti\  =  Xa [n  +  1]  —  Xa [n\ 

finance,  turbulence,. . . : 

Curvature (("(g))  =  A2  <<  1 

If  A2  <<  1,  XaW  “is  close  to”  rand,  walk  Ra[^] 

SRa[ti]  =  e[n]eA^A^/2, 

n\  :  white  noise,  DaM  :  known  Gaussian  process. 


Theorem  9  (Meaning  of  “is  close  to”) 


2  In  M([(n  —  1)A,  nA])  1  law 


A 


n 


A^0+ 


{0A[n]j 


n 


Moreover, 


•  E  (In  \6XA[n]\q  . . .)  =  E  (In \5RA [n] \q  . . .)  (1  +  o(A2~e)) 
and  when  E  (5X A[n]g  . . .)  <  -foo 


E  (SXA [n]q  . . .)  =  E  ( 5RA[n]q  . . .)  (1  +  o(A2-£)); 
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Application  to  Parameter 
Estimation 

(Bacry,  Kozhemyak ,  Muzy ,  2004) 


Cov(\n\SXA[0}\M\SXA[r}\) 
zz  Cov (In  \5Ra [0]  | ,  In  |5i?Ab~]|) 

—  Var (In  | e|)  +  \2Cov(Qa[Q],  Oa[t])/4 
~  —A2  In (t/T),  for  1  <  r  <  T 
=^>  independent  of  A 

Estimation  of  the  multifractal  coefficient  A2 
— »  linear  fit 

Integral  scale  T 
— »  decorrelation  scale 

Variance  a2 

— *  variance  of  the  increments 


GMM  estimation 


30 


Parameter  Estimation 
(changing  A) 

(Bacry,  Kozhemyak,  Muzy  2004) 


S&P  500  intraday  data  (5mn  ticks,  1996-1998) 


:  5mn  returns 
:  30mn  returns 
:  lh  returns 
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Other  applications 

( Bacry ,  Kozhemyak ,  Muzy  2004) 


SXA 

n 

~  6 

n 

^AOaN/2 

Variance  prediction/estimation  (Wiener  filtering) 
•  Prediction  of  fi[n] 

QaM  =  /12  *  hi  (|£Va|)  [n]  (J12  causal) 

— »  MLE  of  exn*W 


Value  at  Risk  prediction/estimation 
•  estimation  of  v(p)  such  that 


Prob  {\5X/s[n\  \  >  v(p)  \  \  SX^[n  —  fc],  k  >  0}  =  p 


Edgeworth  expansion 
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